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Abstract. We consider the problem of picking a physically motivated vacuum 
state on a spherically symmetric spacetime with an extra conformal Killing vector, 
as opposed to an extra Killing vector as in the Schwarzschild case. Considering a 
conformal symmetry instead of a symmetry allows us to consider spacetimes that 
are dynamical and not static (like Schwarzschild). The extra conformal symmetry 
allows us to calculate the response of particle detectors however. We look at the 
specific example of a self-similar LTB spacetime that represents a spherically 
symmetric but inhomogeneous cosmology. We remark that the above procedure 
might be applied to a spherically symmetric collapse solution that represents 
black hole formation so that one can calculate the detailed spectrum of Hawking 
radiation during a collapse. 
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1. Introduction 

In order to calculate the response of particle detectors following a given trajectory in 
a dynamic spacetime, one often picks spacetimes that have a large amount of sym- 
metry. Symmetries of the spacetime make PDE's such as the Klein-Gordon equation 
separable, allowing one to give a complete set of solutions to the PDE in terms of solu- 
tions to several ODE's. The correlation function can be obtained from the complete set 
of solutions and hence the response of an Unruh-Dewitt detector p] can be calculated. 

The response of detectors in anisotropic Bianchi spacetimes can be calculated 
due to the existence of three translational Killing vector fields [2j. In the eternal 
Schwarschild spacetime, existence of spherical symmetry and a global time transla- 
tion Killing vector field allow these calculations to be done everywhere [2J, [3]. 

In [3] black holes without the global time translation Killing vector field were 
considered. It was demonstrated that all black holes formed from collapsing matter 
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radiate with the characteristic black body spectrum. However, the trade off for mak- 
ing an argument for black hole solutions formed from a general spherically symmetric 
matter is that the spectrum of radiation can only be calculated on future null infinity. 

The problem with requiring more symmetry is that it often makes the physical 
system less realistic. For example, the global time translation Killing vector in the 
Schwarzschild spacetime specifics a black hole that was not formed at a finite time 
in the past. Hence, by assuming this symmetry one would not be able to obtain a 
physically realistic response of particle detectors away from future null infinity for 
black holes formed from collapsing matter. In the cosmological realm, homogeneity 
is assumed to be true only at the largest scales and it is likely that at smaller scales 
the geometry of the universe is more like that of a Swiss Cheese model [5]. The 
swiss cheese model contains inhomogeneous regions such as a Lemaitrc-Tolman-Bondi 
(LTB) spacetime. 

We give a physically well motivated prescription for choosing a vacuum state 
that can be applied to any conformally coupled quantum field in a general spheri- 
cally symmetric spacetime with an extra conformal Killing vector field. The reason 
for considering this type of spacetime is that it is more physical than considering a 
spherically symmetric solution with a time-like Killing vector (Schwarzschild) which 
is completely non-dynamical. The late time response of comoving detectors is given 
for the self-similar LTB spacetime given in [5] . 



2. A self-similar LTB spacetime 

This section is an overview of the results given in [5] that are needed for this paper 
(for a full discussion one should consult [5]). 

A model for an inhomogeneous, spherically symmetric universe filled with 
pressure-less dust is given by the following LTB metric. 

ds 2 = dt 2 - cxp(2A)dr 2 - r 2 S 2 dn 2 (1) 

where A = A(i, r) and S = S(t, r). The stress energy tensor is 

Tab = pu a U b (2) 

where p is the density and u a is the 4- velocity of the dust . 

To solve the Einstein equations a self-similarity is imposed by making A = \{t/r) 
and S = S(t/r). The self-similarity implies the existence of a vector field (£ n ) satisfying 
the equation 

£a||6 + £&||a = 2 .9o6 (3) 

where ||6 denotes covariant differentiation. £ a is called a homothetic Killing vector 
field. Self-similarity and the Einstein field equations imply that S and A satisfy the 
following equations 

S' 2 =2E + 2/S (4) 
exp(A) = (S - sS')/VT+2E (5) 
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where s = tjr and S' — ^j. E is an integration constant representing the energy of 
the shell in the infinite past. The case E = will be used. The solution to equation 
(JD) in the E = case is 

s + a = ±^S 3 / 2 (6) 

where a is an integration constant representing the degree of inhomogeneity in the 
metric. 



In the t, r coordinates the line element appears to be singular at r = [5]. 
Since the behaviour of the Klein-Gordon field will be needed near r = 0a coordinate 
transformation will be made to the coordinate to given by r = 2w 3 /2ip, where to is 
some arbitrary time scale. In the tu coordinate the line element is 

ds z = dt 2 



to 



t + 3-juj 



t + 7CJ 3 



du z + oj z dW 



(7) 



where 7 = 6a/(27i§) 



It is generally true that for a spacetime with a conformal vector field (CKV) there 
exists a conformally transformed spacetime in which that C.K.V. becomes a Killing 
vector field. To prove this consider a spacetime metric with a conformal Killing 
vector field satisfying 

Lk9ab = (j>{x c )g ab (8) 

where is the Lie derivative along the trajectory of the vector field k a . Consider a 
parameter r along the trajectory of the vector field k a so that k a can be represented 
by the derivative -4^ . The conformally transformed spacetime given by 



jab = cxp ^- J dr0(r)^ g a , 



g a b = cxp ( - / dr0(r) I g a b (9) 



is a solution of the equation 

Lk9ab = (10) 
confirming that k a is a Killing vector field in the spacetime with metric g ab - 

A coordinate can be chosen along the trajectories of the Killing vector field k a 
(say r) so that the metric is independent of this coordinate. In a spherically sym- 
metric, 4 dimensional spacetime this guarantees that the Klein-Gordon equation is 
separable and hence one can solve for a complete set of solutions. 

The spacetime described above has a conformal Killing vector field given by 

C d^" t dt + 3d^ { ' 

The coordinate y will be chosen so that ^ is the conformal Killing vector field. The 
relationship to the t and co coordinates can be obtained by equating (jllj) and the 
relation 

d dej d ^ dt d 
dy dy 8uj dy dt 
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We obtain 

w = Acxp(y/3) t = Bcxp(y) (13) 

where A and B are integration constants that are independent of y. We choose A = to 
and B = t^x, where x will be another new coordinate. The coordinate x is constant 
along curves generated by the self-similarity transformation. 



The line element in the y, x coordinates is given by 



ds 2 = tl exp(2y) 



dy 2 -JL( dx+ * d y) 2 + xV/W 
a z/i 3 



(14) 



where a = 1 + 7a; 3 and 6=1 + 37a; 3 with 7 = 6a/27. The exp(2y) factor is the only 
dependence of the metric on y. Therefore, d/dy is a conformal Killing vector in the 
above metric g ao and a Killing vector in the metric g^ given by 

g a b = cxp(~2y)g ab (15) 

The conformal Killing vector d/dy becomes null (g(d/dy,d/dy) = 0) when the 
coefficient of the dy 2 term in ([14")) vanishes. The surface where the C.K.V. d/dy 
becomes null defines a conformal Killing horizon [6]. So there is a conformal Killing 
horizon at the solutions of the following equation 

In order to count the number of positive roots in x of the above equation we first 
factor the left hand side using the difference of squares. 

b x\ ( b x\ , 

The factor on the right is always positive for 7 > and therefore has no roots for 
x > 0. Making the substitution z = x 3 into the factor on the left and simplifying gives 
that the positive roots of (|T5|) has the same positive roots as the following equation. 

27 7 3 2 4 + 277V + 9 7 z 2 + (1 - 27-/)z - 27 = (18) 

There is only one sign change in the coefficients of the polynomial on the left hand 
side of (fTSl) . Therefore there can be at most one positive root. Using z — x 3 we see 
that there is only one positive real root for x > and hence only one conformal Killing 
horizon in the expanding, inhomogeneous universe (t > 0). 



We now introduce the mathematical tools to study the evolution of a bundle of 
null geodesic curves in the above spacetime. The purpose is to find other types of 
horizons. The scalar expansion of a congruence of geodesic curves is defined to be 

= «fi„ (19) 

where u a is a tangent vector to the congruence of curves given by d/dX where A is an 
affinc parameter along the curves and | \a is the covariant derivative with respect to the 
coordinate x a [7] . This quantity defines how the cross-sectional area of a small bundle 
of curves changes. For example, a bundle of curves representing radially outgoing 
light rays in Minkowski spacetime would have positive expansion because the curves 
move away from each other as they move to increasing radius, while the expansion for 
radially ingoing light rays in Minkowski spacetime would be negative. 
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A trapped surface of a spherically symmetric spacetime is defined as the surface 
where the expansion is negative for both ingoing and outgoing sets of null geodesies 
[5]. The apparent horizon is the boundary of all the trapped surfaces, meaning that 
the expansion is zero for one set of radial null geodesies and negative for the other. 
For example, the event horizon of the Schwarzschild spacetime is an apparent horizon 

i- 

It can be shown that the conformal Killing horizon (CKH) is also an apparent 
horizon of the spacetime (0 by calculating the expansion of the congruence of the 
two independent, future-directed, radially moving, null geodesies. The ingoing and 
outgoing tangent vectors, denoted by k a _ and fc+ respectfully, are calculated by solving 
the equations 

k b k? {b = (20) 

k a k a = (21) 

For convenience we will work in the conformally transformed spacetime g a \, given by 
(|15|) . We can work in the conformally transformed spacetime because null geodesies 
remain null geodesies under conformal transformations. The tangent vectors in 
afhnc paramctrization in the original spacetime can be obtained by a conformal 
transformation even though working in the barred spacetime (g a b) changes the 
geodesies from afhne paramctrization to non-affine paramctrization. Taking the null 
geodesies to be defined by an affinc parameter A in g a b and working in the coordinates 
of ([7]) we obtain the following expression for (|2"Tj) 

y-^ (i+ (* + ^ (* + !*)) = c 22 ) 

with 

fc a = (y,±,0,0) (23) 

where y = dy/dX. Setting each of the factors to zero gives us an equation representing 
one of the two independent null geodesies. 

y±jL(i + !#=() (24) 

The geodesic equations for the null vector ([2U]) reduce to the following expression 

y-^ ±+ ^ ) = c (25) 

where C is an integration constant. Solving equations (|24|) and (|25[) we obtain an 
expression for in the conformally transformed spacetime 

/ a i/3 \ 

k a ± = C f {l + xb/Za}/3)-\ T — ,0,0) (26) 

It can be verified that in the original spacetime the tangent vectors given by the 
derivative with respect to an afhne parameter are given by 

,1/3 

,0) (27) 



k a ± = Ce- 2y ({l + xb/Sa 1 ^)- 1 ,^^-^ 



It must also be determined where in the spacetime that the null vectors k a are 
future pointing: the reason is because if we find that the expansion associated with 
fc° is negative for example but it is past pointing then as time runs forward the 
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expansion is actually positive. To accomplish this the fc° vectors are expressed in 
the (t, uj) coordinates of ([7]). This is done because both the vectors d/dy and d/dx 
coordinates are time- like in certain regions while d/dio is space-like everywhere and 
d/dt is time-like everywhere. The d/dt component of the null vectors are 

k ± = ct T±xl^ (28) 

So while the component is always positive for x > 0, the component kP_ becomes 
negative when 1 — xb/Za 1 ^ < 0. Hence the vector ki. is always future pointing but 
k°L becomes past pointing when 1 — xb/Sa 1 ^ = 0. Note that this is also the value of 
x where the CKH is found. 

The expansion of each set of null geodesies is found to be given by the expressions 

9 ± = 2e2 ^ T3al/3 ) (29) 
± (3a 1 / 3 + xb)xa 2 / 3 { ' 

It can be verified that for 7 > 8/27 and x > 0, 9 + < and 6>_ > everywhere. 
However k°L becomes past directed at the CKH so the true expansion associated with 
these null geodesic congruences becomes negative. Hence, the CKH is also an apparent 
horizon for 7 > 8/27. This will be of importance later when we use the existence of 
the CKH to give a low frequency cut-off for the two point function. 

3. Conformally coupled Klein- Gordon field 

A conformally coupled Klein-Gordon (KG) field in 4 dimensional spacetime satisfies 
the equation 

d a {V=gg ab d b ^) + \rtI> = (30) 



where R is the Ricci scalar given by 

4 27 7 2 x 6 + I87X 3 + 1 
^="3 ab (31) 



The conformally coupled KG field has the property that if ip satisfies (|30|) in the 
metric g a t then the field ip = (f)~ 1 ip will satisfy the KG equation in the conformally 
transformed metric g a b = 4> 2 gab- We can solve for the conformal fields in the metric 
gab given by (|15j) because the KG equation is separable there and then obtain the 
fields in the metric g a b by a conformal transformation. 

In the metric where d/dy is a Killing vector, the following ansatz is used. 

i>uim = exp(ivy)f vt (x)Y lm (6, cp) (32) 

The vacuum of the ansatz (|32[) is invariant under translations along the Killing vector 
d/dy (because the modes are multiplied by a phase under shifts in y). Substituting 
(|32[) in the KG equation the following differential equation for f v i follows. 

^*.«- + f»W4ff' 'fW^ ,33, 



b 2 J dx 2 \ 3 dx \ \ 9 b 2 J J J dx 

-v 2 abx 2 + J-^(abx 3 ) + 1(1 + 1)-^ - ^ 2 (27 7 V + I872; 3 + 
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The singularities at x = are regular; the powers of x appearing the equation 1331 
have a finite negative power. Therefore (|53"|) can be solved perturbatively by expressing 
the coefficients as a series expanded around x = 0. The two independent solutions 
take the form 

oo 

Ui = * l Y. anx " ( 34 ) 



oo 

fui=x- l - 1 J2b n x n (35) 



where I > 0. The solution (|34[) is regular at x = while (|35]l is singular there. The line 
clement (|T4")) at x = reveals that the metric is regular there. The d£l 2 term in ((H)) 
is multiplied by x 2 , so it docs not contribute at x = 0. The same behaviour occurs 
in flat spacetime in spherical coordinates at r = even though it is a regular point, 
therefore we expect the spacetime to be completely well behaved at x = 0. Since the 
metric is well behaved at x = 0, it is expected that the modes will be also. Hence, 
only equation (f3~4")l will be used for the f v i part of the modes. 



Requiring a vacuum that is invariant along the Killing trajectory reduced the 
arbitrariness down to two complex constants (bo and ao in (|34j) and (|35jl ). Imposing 
regularity at x = implied bo = 0, reducing the choice of vacuum down to one complex 
constant (ao). Finally, the canonical commutation relations of the KG field operator 
imply the modes must satisfy the following normalization. 



(36) 



(VVZ'm',VVZm) = / dx \f ZI gg a (ipt'l'm'dalpvlm ~ ^vlmda^t'V m') 

= i5(v - v')5u>6 mm < 

where the integration is done over the x 1 , x 2 and a; 3 direction. The normalization 



(|36| fixes one of the two remaining real constants leaving only a phase factor unfixed. 
The phase factor disappears when calculating correlation functions. Hence imposing 
the conditions of regularity at x = 0, invariancc under y translations and the normal- 
ization (|36p uniquely specifies the vacuum state. 

Using the ansatz (|32p the normalization condition (|3"6"|) reduces to the following 
condition on ao- 

-l 



| a.o 1 2 = i ( 2ftabx 2 
where 



2iujx 



^3 

b 2 



)( 



h*,ih' u i — h u ih*,i 



Ki = x l 1 1 + ]T 



(37) 



(38) 



The integration (|3^|) is done over a constant x surface so the condition (|3T|) is evalu- 
ated at some constant x = u. 



The correlation function in the vacuum state defined above (the vacuum state 
which is regular at r = and has a conformal symmetry along d/dy) can be calculated 
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poo 

G(x,x') = y2 / dojcxp(- i Lo(y-y')-(y + y'))f^(x)f: i ( X ')Y lm (e^)YC m (9 , , V ' 



)(39) 



Only the correlation function evaluated along the points 9 = 9' and tp = tp 1 will be of 
interest. The following formula can be used to simplify the sum over modes: 

m—l 

Y lm (9^)Y l * m (9^) = (2l + l)/4n (40) 



m=—l 

giving 



OO POO r>/ _1_ 1 

G(x,x')=J2 faew(-Mv-v')-(v + l/))M*)f:i&)-7— (41) 

1 n •* 



1=0 

|2 



Evaluating |do| f° r a ^ the lowest order in x, the following expression is obtained 

2 ^ 
= 47r(l+ £ ,u 3 )(l + 3 5 M 3 )ii 3 a; ^ 

where u is the constant x surface that the modes are normalized on. The Taylor 
expansion for f w i is given by 

f ul = a x l (l + O(x 2 )) (43) 

It can be seen that the correlation function will be divergent at u = due to the 
normalization on the modes. To solve this an IR cutoff must be imposed. A value 
for the cutoff is given by the constant x surface that is the conformal Killing horizon. 
This surface is also an apparent horizon, so it would be natural to impose that no wave 
mode have a wavelength that is larger than the apparent horizon. The single solution 
to (|16p is denoted by x = j3. The function f w i{x) has approximately the wavelength 
2tt/u>. Therefore only modes with w > 2tt/[3 will be included in the sum to calculate 
the correlation function. 



(44) 



f°° 21 + 1 
G(x,x') = 2_. / dujcxp(-iuj(y- y') - (y + y'))f u l(x)f*i(x') — ■ 

Performing the integration and summation gives 

G(x, x') = N(u, 7 )^ ^ (45) 

exp(y + y') 

where N(u, 7) is only dependent on the normalization surface x = u and the 
inhomogeneity parameter 7. Ei(l,x) is the exponential integral special function and 
P(a, b, x) is the associated Legendre polynomial of the first kind. They are given by 
the following. 

/•oo 

Ei(l,x) = dy exp(-yx)y" 1 (46) 



TOW )- 1" 1 ™?^:™- 1 ™ («) 
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The correlation function is being considered for small values of x so it will be 
expanded around x = in a series. Keeping only the leading term the following is 
obtained 



The above function is the correlation function for an observer following the trajectory 
x = 0. It is also the leading order approximation to the correlation function for an 
observer at small x. 

One might expect that the modes that define the vacuum state should oscillate 
infinitely rapidly as the apparent horizon is approached. For example, the modes 
defining the Unruh vacuum state oscillate infinitely rapidly near the future horizon of 
the eternal Schwarzschild spacetime. This is not the case for the other two standard 
vacuums on the Schwarzschild spacetime. However the Unruh vacuum is thought to 
be the physically correct vacuum for an object collapsing to form a black hole. In 
other words, it is the correct vacuum for a star that collapses to form a future horizon. 

Solving (|33|) numerically for regular initial conditions at x = 0, it is found that 
the modes fu>i{x) seem to oscillate infinitely rapidly at the conformal Killing horizon. 
Hence the modes ([5!?]) that specify the vacuum state have the physically desired 
property that they oscillate infinitely rapidly at the conformal Killing horizon. 

4. Response of Particle Detectors 

The response rate of an Unruh-Dewitt detector will be calculated for observers co- 
moving with the dust of the above LTB metric. The expression for the response rate 
is given by 



s is the amount of proper time passed since the initial time tq at which the detector 
is in its ground state. G(x, x') is the correlation function and r is the time at which 
the response rate of the detector is evaluated. 

The r and t coordinates of the metric ([1]) will be used to parametrize a co-moving 
observer, r is constant for a radially infalling dust particle and t is the proper time of 
this dust particle. The coordinate transformation {x, y) — > (t, r) is given by 



The trajectory can therefore be parametrized by the y coordinate instead of t by 
using the relationship y = exp(t). Making the substitution y = cxp(r — s) the integral 
becomes 

F v (fj) = dy' eM-^Wr^Eiil, j(e + i(y - y'))) exp(-(y + j/))) (51) 

where yo = exp(ro) and y = exp(r). The correlation function for the observer at x = 
given by (|48p was used. 




(48) 




(49) 




(50) 
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The response rate of a particle detector near x — approaches zero as the proper 
time of the detector advances. The null response of the detector as y —> oo can be 
demonstrated by noting that the divergent part of the correlation function is due to 
the small x behaviour of the function Ei(l, x), that is: 

Ei(l, x) = -7 - ln(x) + 0(x) (52) 

Using the expansion (|52|l in ([5"Tj) it can be seen that the integral is finite for fixed yo 
and y. As y — > oo the integral (|S"Tj) (without the exp(— y)) factor is finite. Including 
the exp(— y) factor, it is easily seen that the integral is zero. Hence the detector does 
not click as it moves forward in time away from a fixed yo- 

A similar argument can demonstrate that any comoving detector (constant r, 9, <f> 

1/3 

trajectory) will not click. This can be seen by using the relation x = (9r/2t) ' . As 
t — > oo, x — > 0. The integral ([5*Tj) can be divided into two parts, one for late times or 
x < e, and on for early times x > e. The early time part is finite and the late time 
part is finite as t — > oo as demonstrated above. Hence, both parts are exponentially 
damped by the conformal factor exp(y) as in the small x case, and the detector does 
not click at late times for any comoving observer. 



5. Discussion 

The selection of a vacuum state is a non-trivial process for a general background met- 
ric. In the eternal Schwarzschild spacetime, there are three common choices for the 
vacuum state. The Hawking-Hartle vacuum [9] is a state that is regular on the past 
and future horizons. The Boulware [TU] vacuum is symmetric with respect to the time 
translation Killing vector field d/dt. Finally the Unruh vacuum state [1] is symmetric 
with respect to the Killing vector on the past horizon and symmetric with respect to 
the Killing vector on past null infinity. 

We have given a physically motivated prescription for choosing a vacuum state 
of a conformally coupled quantum field on the spacetime given in [5]. The response 
of a particle detector was given in this state so that the response along any trajectory 
could be given by integrating numerically. The late time response of comoving particle 
detectors was found to be zero in this vacuum. 

The above prescription for the vacuum state can be used for any conformally 
coupled field on a spherically symmetric background with an extra conformal Killing 
vector field. In particular, there are solutions to the Einstein Field Equations repre- 
senting spherically symmetric matter collapsing to form a black hole [TTj . The solution 
in [TT] also has a conformal Killing vector field. One could therefore study the detailed 
spectrum of Hawking radiation from this black hole. 

A solution of Einstein's equations representing matter collapsing to some central 
point (r = 0) should be regular along the line of points specified by r = before 
some specified time when the singularity forms. The metric is regular at these points; 
hence we expect the vacuum state to also be regular there which gives us one of the 
conditions we have used to specify the above vacuum state. Further, we expect a 
conformally coupled field to have the same conformal symmetry as the background 
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spacetime, which gives us the second condition on our vacuum state. These two con- 
ditions uniquely specify the vacuum state and for the above reasons we believe that 
it is the physically correct vacuum for spherically symmetric collapse solutions with a 
conformal Killing vector. 

We note that the condition of regularity at r = cannot be used on the eternal 
Schwarzschild spacetime because there is always a singularity in the metric at r = 
and therefore there is no reason to expect that other scalars (such as a Klein-Gordon 
field) would not also blow up there. By 'eternal Schwarzschild spacetime' we mean 
the spherically symmetric spacetime that everywhere has a time translational Killing 
vector. This spacetime cannot represent matter collapsing to form a black hole. The 
conditions we have used to define a vacuum are more physically realistic in the sense 
that they apply to a spacetime representing a black hole that was formed at some 
finite time in the past. 
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